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Abstract 

We report observations of the lunar helium exosphere made between Decem- 
ber 29, 2011, and January 26, 2012, with the Lyman Alpha Mapping Project 
(LAMP) ultraviolet spectrograph on NASA's Lunar Reconnaissance Orbiter Mis- 
sion (LRO). The observations were made of resonantly scattered He I from 
illuminated atmosphere against the dark lunar surface on the dawn side of the 
terminator. We find no or little variation of the derived surface He density with 
latitude but day-to-day variations that likely reflect variations in the solar wind al- 
pha flux. The 5-day passage of the Moon through the Earth's magnetotail results 
in a factor of two decrease in surface density, which is well explained by model 
simulations. 
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1. Introduction 

We have recently reported on the detection of helium in the lunar exosphere 
using the Lyman Alpha Mapping Project (LAMP) ultraviolet spectrograph on 
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NASA's Lunar Reconnaissance Orbiter Mission (LRO) (IStern et all 120 121) . These 
were the first measurements of He since the in situ Lunar Atmosphere Composi- 
tion Exp eriment ( LACE) measurements made follow ing the Apollo 17 mission 



in 1972 (|Hodgesl 1 19731 : iHodges and Hoffmanl 11974 ) and the first time that He 



has been seen by remote sensing instrumentation. LAMP measured resonantly 
scattered solar He I X584 from above the sunlit lunar limb. To correct for the in- 
terplanetary He I A.584, the spacecraft executed a slew to point LAMP at fixed 
celestial coordinates keeping the sky background constant during the maneuver. 
These observations provided a snapshot of the He distribution at just a few specific 
times, and viewing across the poles. 

On December 14, 2011, the LRO Project altered the spacecraft's polar orbit 
from a near-circular 50 km altitude to a more stable, fuel efficient elliptical orbit 
with apoapsis near 200 km and periapsis at 30 km. This change occurred when 
the orbit was very close to the terminator (at the terminator the [3-angle, defined 
as the difference between the Sun direction and the spacecraft orbit plane, is 90°). 
As soon as [3 decreased to the point where it was feasible to observe the dark 
lunar surface with the full LAMP aperture, we found that the enhanced He col- 
umn along the line-of-sight from apoapsis (on the night side the He scale height 
is about 130 km) made detection of He I X584 fairly straightforward because of 
the absence of sky background. Moreover, it would allow continuous monitoring 
of the He atmosphere until (3 reached the point where the spacecraft was totally 
in shadow, or for roughly a one-month period. The technique of observing illu- 
minated atm osphere above the dark lunar surface is the same as that employed by 



Fastie et al.l(|1973|) with the Apollo 1 7 Ultraviolet Spectrometer, except that in that 
case the near-equatorial orbit permitted only a few minutes of observation at each 
terminator crossing. 

In this paper we present the details of the spectroscopic analysis of the ad- 
ditional, more extensive observations made in late December 2011 and January 
2012, discuss the calibration issues, and describe the parameterization used to in- 
fer the He surface density from the data. We will also discuss our results in terms 
of the solar wind alpha particle flux and recent modeling efforts. 

2. Observations 

LAMP is a lightweight, low-power, imaging spectrograph optimized for far- 
ultraviolet (FUV) spect roscopy of the nightside lunar surface and atmosphere 



(|Gladstone et al.U201 0b). It is designed to obtain spatially resolved spectra in the 



570-1960 A spectral band with a spectral resolution of 34 A for extended sources 
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that fill its field-of-view. The slit is 6.0° long, with a width of 0.3°. Each spatial 
pixel along the slit is 0.3° long. From an altitude of 50 km the spatial resolution 
is 0.26 x 0.26 km 2 . The use of LAMP for the detection of gaseous emission was 
demonstrated by the observations of the plu me produced by the LCROSS impact 



on October 9, 2009 (Gladstone et all 1201 0a). 



A diagram showing the observing geometry for the "frozen" elliptical orbit 
is shown in Fig. \T\ for two different dates. The initial time in each plot is the 
equator crossing time on the day side. Each orbit, of roughly two hours, generates 
a separate data file of time-tagged (pixel list) photon events. As shown in the 
cartoon of Fig. |2l the orbit crosses the north pole at roughly t+1800 s, descends 
to the night side equator and crosses the south pole near periapsis at ~t+5400 s. 
Near the poles there are bright sunlit spots that cause the detector high voltage to 
be temporarily reduced producing data dropouts, so that our useful data range is 
typically from t+2200 s to t+5100 s, corresponding to a latitude range of 80° N to 
80° S. For a given day, we co-add the spectra from all of the orbits to enhance the 
signal/noise ratio. 

In Fig.[D the shaded area represents the distribution of illuminated atmosphere 
for a nadir-looking line-of-sight. The shadow height, the solid line, is calculated 
for a spherical moon and and is uncertain to the level of a few km. The dashed 
line is the LRO altitude. The figure on the left is for December 29, 2011, with a 
P-angle of 83.2°. The figure on the right is for January 16, 2012, (3=64.7°, and for 
a good portion of the nightime orbit the spacecraft is in shadow. These latter data 
will prove invaluable in the extraction of the He emission signal. 



3. Data Analysis 

3.1. Lunar atmosphere spectrum 

As an example of the short-wavelength region of the LAMP spectrum, Fig. [3] 
shows a composite co-added spectrum from 12 orbits on January 16, 2012, in 
which the count rate in the illuminated rows of the detector are summed for times 
between 2100 and 2600 s after the start of the orbit. This corresponds to the left 
shaded region in the right panel of Fig. \T\ The background is due to grating scat- 
tered Lyman-oc photons, from the interplanetary background reflected from the 
dark surface, and the steep rise to the short-wavelength limit is caused by elec- 
tronic pile-up at the edge of the sensitive area of the detector. Another artifact 
exists at wavelengths 760-800 A. Overplotted in red is a spectrum with the space- 
craft totally in shadow, from 3300 to 4700 seconds. A very small He I A584 signal 
remains, presumably due to reflected interplanetary He I from the dark lunar 



3 



surface. The difference between the two gives the emission from the illuminated 
column, as seen in the lower panel. In extracting the He emission as a function of 
time, the shadow spectrum is used as a template for subtracting the instrumental 
background. 

3.2. Orbit variation 

Fig. \T\ suggests that the observed He I ?i584 emission rate should follow the 
height of the shaded area as a function of time during the night side of a given 
orbit. This behavior is indeed found, as can be seen in Fig. |4] for the two dates of 
Fig. [Q To enhance the signal/noise ratio of these data we have co-added all of the 
orbits for a given date. [A day-by-day compilation of all of the data for the entire 
period is available in the Supplementary Online Material, Figs. [6}{8]] . The data are 
parameterized as described in Section I3~4l 

3.3. Calibration 

Calibration at 584 A is difficult. Although the instrument sensitivity, or effec- 
tive area, was measured pre-flight in the laboratory, there is no convenient way to 
monitor the in-flight performnce at this wavelength as is routinely done for wave- 
lengths longward of the interstellar hydrogen absorption edge at 912 A using a 
selection of standard stars observed on a regular basis. Instead, we chose to rely 
on an observation of strong interplanetary He I X584 emission, made on July 20, 
2011, when LAMP was viewing towards the dark limb with the Sun only several 
degrees below the limb. The line-of-sight viewed a bright region of the interplan- 
etary He focussing cone. We assessed this background and its variation on the 
sky using a standard helium "hot" model ( Ajellol Il978|) that we have tuned by 



fitting unpublished Galileo EUVS interplanetary cruise observations of interplan- 
etary heliu m obtained at the s ame time as published observations of interplanetary 



hydrogen (|Pryor et al.L |1996|) . Model parameters used for the interstellar helium 
at large distances from the Sun are a density of 0.015 cm~ 3 , velocity of 26.3 
km s , and a temperature of 6300 K. The solar He I X584 li newidth is taken 



as 13 mA full- width at half-maximum (Doppler width 78 mA, Lallement et al. . 



20041) . Solar He I X584 brightness values are taken from the Thermosphere Iono- 



sphere Mesosphere Energetics and Dynamics (TIMED) satellite' s Solar EUV Ex 



periment (SEE) Level 3 data available online (|Woods et al.L 1 19941) . He atom pho 



toionization rates are taken fro m the Space Environ ment Technologies online So 



lar Irradiance Platform v2.37 (|Tobiska et all |2000J). We integrated over the slit 



position on the sky over the five minute duration of the observation and deduced a 
sensitivity at 584 A of 0.75 counts s~ ! rayleigh -1 . We have verified consistency 
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by comparing nighttime observations of the sky with the models for those partic- 
ular dates and view directions. All of the nadir observations have emission rates 
of the order of or less than 1 rayleigh. 

3.4. Data Modeling 

We parameterize the He atmosphere a s a surface bounded exosph ere using the 
classical model of IChamberlainl (|1963|) . lHartle and Thomasl (|1974|) have shown 
this to be a valid approximation for He at low altitudes, such as those considered 
here. We tak e the exobase tempera ture to be the surface temperature which we 
set at 100 K (|Vasavada et all |2012|) . leaving the helium density at the surface as 
the sole variable. We assume that He I X584 is optically thin and calculate a flu- 
orescence efficiency (g-factor) using daily solar He I X584 fluxes from SDO/EVE 
dWoods et all |2010() and the solar linewidth quoted above. The density model is 



then integrated between the spacecraft orbit and the shadow height to calculate a 
predicted brightness which is then fit to the data using a linear least-squares algo- 
rithm, to infer the He surface density. The fits are illustrated in Fig. |4] and in the 
Supplementary Online Material Figs. (ME 

For all of our dates, the fit of the model to the data is excellent, and implies 
that the He surface density is roughly constant with latitude. Because of the night 
time He scale height of ~ 130 km near the dawn terminator, there is also not much 
sensitivity to the choice of surface temperature, as a variation of ±20 K gives only 
a ±9% and ±2% change in calculated brightness from LRO altitudes of 200 and 
35 km, respectively, which is within the range of statistical uncertainty in the data 
fits. 

3.5. Temporal Variability 

The variation of the inferred surface He den sity w i th date is shown in Fig. [5J 
During this period, the longitude, defined by iHodgesI (| 19731) as the angular dis- 
tance from the sub-solar meridian, is decreasing. At the dawn terminator, the 
longitude is 180°+|(3|, and this is shown on the figure. This quantity is analo- 
gous to local time and is useful in comparing the data with models in the liter- 
ature. The mean density inferred for longitudes within 15° of the terminator is 
~ 2.0 x 10 4 cm~ 3 , in good agreement with the Apollo 17 LACE data taken at a 
latitude of 20° for the same range of longitudes. On a given date, orbit-to-orbit 
variations are less than the statistical uncertainties in the count rates, so that the 
time scale for density variability is greater than 2 hours (the orbit period), but of 
the order of a day, the time scale for thermal escape of He ( Hodgesl Il973|) . Our 
data show a variation of a factor of three in derived density over a 28 day period. 
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One striking feature stands out: a decrease in the surface He density by a factor 
of two in a five-day period beginning on January 7, 2012. This decrease occurs 
during a period in which the Moon passes through the Earth's magnetot ail. The 
shape and magnitude of the decrease match very well the prediction of iHodges 
(I1978|) . and are discussed below. 



4. Discussion 



Following the Apollo 17 mission, discussion and modeling of the LACE mea- 
surements (that were made over 10 consecutive lunations) focussed on the ques- 
tion of the efficacy and distribution of solar wind alpha particle neutral ization, 



ther malization, and eva poration as th e primary source of the He exosphere (IHodges and HoffmanL 
Il974|) . Modeling by iHodgesI Jl973 ) was able to reproduce the observed lon- 
gitudinal dependence and magnitude of the LACE measurements, with the as- 
sumption that these processe s were almost 100% efficient. Recent modeling, by 



Leblanc and Chaufrayl (1201 1|) . supports these conclusions. The temporal variabil- 



ity observed by LAMP would then be expected to correlate with the solar wind 
alpha flux. Similar variations in the night time He density were seen in the LACE 
data and shown to correlate with th e geomagnetic Kp index, a gain suggesting a 
dependence on the solar wind flux (|Hodges and Hoffmanlll9 74). 



Using the Monte Carlo lunar exosphere model of [Killen et al.l (I2012|) . we sim- 
ulate the time-dependent helium density in the lunar exosphere. We assume a 
source of solar wind alpha particles that become neutralized and thermalized 
through interaction with the lunar regolith on short time scales. The flux incident 
on the lunar surface falls off as a function of solar zenith angle. The aberration 
in the apparent solar wind direction of 5° is neglected. The neutral helium atoms 
are released into t he atmosphere with a velocity taken from a Maxwell-Boltzmann 



flux distribution (Smith et al 



Diviner data ( Vasavada et al 



19781) at the local surface temperature taken from 
20121) . The model tracks the helium atoms as a 



function of time assuming that the outbound velocity comes from the local ther- 
mal distribution after each encounter with the surface. Tracing the evolution of 
the helium exosphere with time enables a time-dependent source rate to be imple- 
mented. 

Upstream solar wind monitors can provide measures of the solar wind speed 
and the alpha particle density at the LI point. However, until calibrated data are 
available for this time period, we assume a constant solar wind alpha flux of 1 .2 x 
10 7 cm~ 2 s _1 , except when the Moon passes through the Earth's magnetotail. 
We do not use a different alpha flux for the magnetosheath than for the unshocked 
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solar wind. Using data from the twin Acceleration, Reconnection, Turbulence, and 
Electrodynamics of the Moon 's Intera c tion w ith the Sun (ARTEMIS) spacecraft in 
orbit about the Moon (lHalekas et al. . 2011 ). we identify the times around full 
moon when the ion flux at the Moon is drastically reduced compared to nominal 
solar wind. During these times, we set the source rate to zero. 

In Figure [5j we show the model density (red diamonds) as a function of time 
compared to the LAMP observations (black asterisks). The model density was 
taken at the lunar surface from the longitude of LRO's nightside passage for each 
day. The trend of increasing density at the beginning of January in the model is the 
result of the progression of LRO's orbit further from the terminator. Deviations 
between the model density and the LAMP observations are probably due to short- 
term variations in the solar wind alpha flux. The decrease in He density beginning 
on January 7, however, does not appear to correlate with a comparable decrease 
in solar wind flux. As noted above, this time corresponds to the Moon's path 
entering into Earth's magnetotail which shields the surface of the Moon from 
the outflowing solar wind. During this time the He density appears to decrease 
exponentially with a time constant of ~5 days which we associate with the thermal 
loss rate of He atoms at the corresponding surface temperature. Once the Moon 
exits from the magnetotail, the density returns to its former value with a similar 
time constant. Our model closely tracks the temporal behavior of the inferred He 
surface density and is also in qualitative agreement with the initial predictions of 



Hodges! (|1978|) . The LAMP data are the first such measurements to confirm these 



predictions. 

We note that a large solar coronal mass ejection (CME) occurred on January 
23-24 and produced an extremely large increase in solar wind flux but also gen- 
erated a large flux of MeV electrons that produced a highly elevated background 
count rate on the LAMP detector and masked any He I X584 signal. After the 
event, the inferred He density does not appear higher than before the event, but 
by this time the (3-angle had decreased to the point where only a small fraction 
of the orbit had illuminated atmosphere and our measurements are limited to high 
latitudes. 

One caveat about our measurement technique is that it does not give us access 
to the complete range of longitudes and latitudes on the night side of the Moon. 
We are limited to range of about 30° in longitude on the night side of both the dusk 
and dawn terminators. However, we expect to be able to measure the dawn/dusk 
asymmetry in He density seen in the LACE data, as well as to track long term 
correlations with the solar wind alpha particle flux. 
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5. Conclusion 

We have described continuous observations of the night side lunar helium ex- 
osphere made over a period of a month in January 2012 using the Lyman Alpha 
Mapping Project ultraviolet spectrograph on NASA's LRO mission. These obser- 
vations show that the inferred surface He density exhibits day-to-day variations 
that likely vary with the solar wind alpha flux and decrease significantly during 
the passing of the Moon through the Earth's magnetotail. Once calibrated solar 
wind alpha particle fluxes measured by the Solar Wind Electron Proton Alpha 
Monitor (SWEPAM) instrument on the Advanced Composition Explorer (ACE) 
satellite are publicly available for January 2012, we intend to carry out a more de- 
tailed investigation of the correlation between the solar wind fluxes and inferred 
surface He density. We also expect additional LAMP data from dawn and dusk 
terminator crossings during the summer of 2012, which we will report in a future 
publication. 
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FIGURE CAPTIONS 



Fig. Q] LRO altitude (dashed line) and calculated shadow height (solid line) for the 
first orbit on December 30, 2011 (left) and January 16, 2012 (right). The origin 
of the time axis is the time of equator crossing on the day side of the orbit. The 
shaded area shows the region of solar illuminated atmosphere below LRO. 

Fig. [21 Cartoon showing the elliptical LRO polar orbit (in red) relative to the 
dawn terminator and the ~ 130 km He scale height. Apoapsis is near the north 
pole (at top), the Sun is to the right. LAMP, viewing towards the nadir, observes 
illuminated atmosphere against the dark surface of the Moon. 

Fig. [3] Average of 12 short wavelength LAMP spectra of the illuminated (black) 
and dark (red) atmosphere on January 16, 2012. The feature centered at 775 A 
is an instrumental artifact. The lower panel is the difference, the average atmo- 
spheric component of He emission. 

Fig. IH Observed count rates for the atmospheric He I ?i584 emission for the two 
dates of Fig. [Q For both dates the data from 12 contiguous orbits have been co- 
added. The error bars are l-o statistical uncertainties. The overplotted models are 
discussed in the text. 

Fig. [5J Surface density of helium is shown as a function of time from LAMP ob- 
servations (black asterisks) and model predictions (red diamonds) averaged over 
each day. Error bars denote estimated 1-a uncertainties in the fit of the data to 
the exospheric column density model. The progression in longitude of LRO's 
nightside groundtrack is provided near the bottom. Times of the Moon's passage 
through the Earth's magnetotail are indicated by dashed lines. 
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Figure 1 : LRO altitude (dashed line) and calculated shadow height (solid line) for the first orbit on 
December 30, 201 1 (left) and January 16, 2012 (right). The origin of the time axis is the time of 
equator crossing on the day side of the orbit. The shaded area shows the region of solar illuminated 
atmosphere below LRO. 



13 




Figure 2: Cartoon showing the elliptical LRO polar orbit (in red) relative to the dawn terminator 
and the ^130 km He scale height. Apoapsis is near the north pole (at top), the Sun is to the right. 
LAMP, viewing towards the nadir, observes illuminated atmosphere against the dark surface of the 
Moon. 
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Figure 3: Average of 12 short wavelength LAMP spectra of the illuminated (black) and dark (red) 
atmosphere on January 16, 2012. The feature centered at 775 A is an instrumental artifact. The 
lower panel is the difference, the average atmospheric component of He emission. 
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Figure 4: Observed count rates for the atmospheric He I A-584 emission for the two dates of Fig.Q] 
For both dates the data from 12 contiguous orbits have been co-added. The error bars are 1-a 
statistical uncertainties. The overplotted models are discussed in the text. 
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Figure 5: Surface density of helium is shown as a function of time from LAMP observations 
(black asterisks) and model predictions (red diamonds) averaged over each day. Error bars denote 
estimated l-o uncertainties in the fit of the data to the exospheric column density model. The 
progression in longitude of LRO's nightside groundtrack is provided near the bottom. Times of 
the Moon's passage through the Earth's magnetotail are indicated by dashed lines. 
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SUPPLEMENTARY ONLINE MATERIAL 



Figures 6-8. The following plots show the same data as in Fig.|4] but day-by-day 
from December 30, 201 1 through January 22, 2012. For each date the number of 
contiguous orbits that have been co-added is indicated in the title, together with 
the day of year (DOY) and year. The error bars are l-o statistical uncertainties. 
The parameters used to fit the overplotted models are also indicated on each plot. 



Fig.|6l December 30, 201 1 to January 6, 2012. 
Fig. [7] January 7, 2012 to January 14, 2012. 
Fig. E January 15, 2012 to January 22, 2012. 
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Figure 6: December 30, 201 1 to January 6, 2012. 
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Figure 7: January 7, 2012 to January 14, 2012. 
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Figure 8: January 15, 2012 to January 22, 2012. 
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